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1. Introduction 

The correct understanding of the Conservation level of a historical masonry monu- 
ment often represents a difficult task to deal with, because of both the high number of un- 
certainties and the difficulties related to the correct representation of the monument’s be- 
havior with simple structural schemes. 

The nowadays museum of Hagia Sophia (which formerly was a church and a 
mosque) is one of the most masonry monuments inherited from past days, thus represent- 
ing an exceptional task for those involved in structural analyses of monumental buildings. 

A correct modeling procedure would require, besides the necessary information on 
the present geometry, the knowledge of all the major Steps in the monument’s history and 
evolution; it is then necessary to acquire information about its original shape and the tech- 
niques used during its erection, about all the events happened during the life of the 
monument (collapses and other significant events), about some possible changes in the 
mechanical behavior (i.e. due to creep or to inelastic deformations). As a matter of facts, 
the present structural behavior is due to a sequence of events which should be known 
when dealing with any intervention aimed at the monument’s protection. 

Some of the aforementioned aspects are enhanced in Hagia Sophia because of the 
peculiar properties of the masonry which has been used during its erection (very thick mor- 
tar beds), because of some important collapses happened during the years, because of a 
series of substantial interventions aimed to the safeguard of the monument (which often 
completely modified its structural behavior). 

It is then clear then the set-up of a numerical model reproducing the behavior of the 
whole monument would be almost meaningless, so that in the present study only a particu- 
lar aspect has been investigated, by using only “partial” structural models. 

In Hagia Sophia, one of the most important aspect is represented by the relevant in- 
clination of both main and buttress piers on northern and Southern sides of the monument, 
as it has been reported in other studies performed by Mainstone that will be described in 
the following [1], The large difference between the values of dorne’ main diameters that 
has been surveyed could then also be due to the “structural weakness” of those structural 
elements; the analyses of displacements occurred at piers’ top level just after their erection 
will then be the aim of the present research. 

After a preliminary section in which some of the most important Steps in the history 
and the evolution of the monument will be sketched, the analysis of piers’ displacements 
will be discussed; the adopted approach consisted in both linear-elastic analyses and 
nonlinear ones, the latter aimed to define “equivalent” linear parameters to be used in 
subsequent finite elements models. 

2. Some historical notes on erection, collapses and rebuilding 

During Nika revolt in January 532 A.D. the existing Hagia Sophia church was fired 
and partially destroyed. Although most of the brick and stone basilica was probably left 
Standing, Justinian decided to create a new cathedral and the first stone was placed just 
40 days after the event. The new church was consecrated on December in the year 537, 



that is less than five years were necessary to build the whole monument, essentiaily the 
same structure as it Stands today. 

For more than nine centuries Hagia Sophia represented the masterpiece and the re- 
ligious center of the Byzantine Empire, the “Megala Ekklesia”; when Constantinople feil to 
the Ottoman army, 28 May 1453, the young Sultan Mehmet II had the Church of Hagia 
Sophia turned into his own mosque. Only in 1934 Aya Sofya mosque was opened as a 
museum and this is its present state. 

2.1 Collapses and rebuilding 

The global aspect of the building at the time when it was consecrated was certainly 
different with respect to the present one. The original dorne probably presented a different 
radius of curvature that was the one of a pendentive dorne, with the same System of large 
base ribs and window openings which is clearly visible nowadays; no one of the buttresses 
and other strengthening works were obviously there. 

As it has been reported on some historical documents, at least five different earth- 
quake events hit Constantinople between 542 and 557 A.D. [2]: on December 14 th , 557, a 
severe earthquake caused serious damages to the dorne so that some restoration works 
were planned. Düring these repairs the eastern semidome, its supporting arch and part of 
the oriental portion of the dorne collapsed, destroying the Lord's Table, the ciborium and 
the ambo (May 7 th , 558). 

A complete re-building of the dorne and of the eastern semidome was decided, and 
internal sides of northern and Southern arches were corrected (by an increasing change of 
thickness) in order to give a more circular shape to the impost of the new dorne. As a mat¬ 
ter of fact North-South and East-West diameters suffered some changes (the former in¬ 
creasing its length and the latter decreasing it) because of the inclination of main piers as 
well as of the thrusts exerted by the dorne. The new dorne was rebuilt steeper and of 
lighter materials. 

Successive collapses in 989 and 1346 were very similar to that of 558, and they af- 
fected structural western elements and oriental ones respectively. After these two last 
events, only damaged parts and close structural elements were repaired. At present time, 
different reconstructions are clearly visible: original portion are the Southern and the north¬ 
ern parts, while western has been rebuild after 989 earthquake and eastern after third par¬ 
tial collapse occurred in 1346. 

The main western arch is thicker with respect to the original one: it withstood the two 
earthquakes in 1343 and 1344 that provoked the collapse of the western one, while only 
small damages occurred during another earthquake inl 894. 

It is likely that, because of collapses and reconstructions, original structures under 
the dorne level were largely deformed maybe even the 557 event but certainly after the 
earthquake in 989. As an example, the four towers in which stairs are located turn out to 
be clearly rotated with respect to main symmetry axes; the rotation, probably due to the 
horizontal thrusts exerted by the main arches, has been hidden by adding external ma- 
sonry layers, as recent restoration works put into evidence. 

2.2 Additions 

Reconstructions altered the original shape of Hagia Sophia, but other three se- 
quences of restoration and strengthening works have been documented through the life of 
the monument: a first phase begun in 1317 in the reign of Andronicus II Palaeologus; a 
second one, under the supervision of architect Sinan, started in 1573; a third phase in 
1847 due to Swiss architects Giuseppe and Giovanni Fossati. It is just worth to remind that 
in 1317 the buttresses at the east wall were erected, while Sinan provided for a series of 



strengthening works and for the erection of a new minaret. Fossati made some changes at 
the dome’s base ievel, introduced some Steel ties and corrected the verticality of some of 
the columns in the exedrae. 

3. Structural Systems in Hagia Sophia 

In Hagia Sophia two main structural Systems can be observed: referring to Mainstone 
book [1] and observing the Van Nice surveying [VN], a primary structural System can be 
identified, which sustains and transfers to lower levels the weight of the superior vaults, 
while the secondary one is composed by inferior vaults and by the galleries. 



Fig. 1 - The Hagia Sophia primary structural System (Mainstone[1]) 
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Both Systems border the main space offered by the central aisle and by the back- 
ground spaces made up by the nartex, the side aisles and the galleries. Starting from the 
upper Ievel of the monument, the structural restraint to the dorne and to the other vaults 
covering the central aisle is essentially due to main and secondary piers, to buttress piers 
(all elements are located outside the central aisle) and to the superior arches connecting 






























































them. The primary structural Systems we are referring to is then composed by these ele- 
ments as well as by the dome and the superior vaults. 

All the other structural elements constitute the secondary structural System, then 
consisting of the arches and the vaults of the nartex, of the major halls offered by lateral 
aisles and of the galleries, besides of their supporting walls and columns not already in- 
serted into the primary System. The secondary System also includes big tympana under 
North and South main arches and walls around externai access ramps. This System is par- 
tially supported by the columns surrounding the central aisle and by the main piers of the 
primary System: a structural dependence of the secondary System from the primary one it 
is then ciear, stating the impossibility of the secondary System of supporting itself without 
the aid of the other structural System. 

The aim of the present work is to focus the structural Situation in 535 A.D., in order to 
understand whether or not the initial configuration of the monument could already be af- 
fected by large dispiacements at the impost level of the dome so that the cross section of 
the dome could not be perfectly circular even at a very early stage. 

In the following, the construction sequence proposed by Mainstone (and nowadays 
accepted by major schoiars) will be adopted in order to give a temporal evolution of the 
building. In the lack of other information, Mainstone follows the historical sequence of 
events as it has been reported by the historian at the Court of Justinian, then reconstruct- 
ing a temporal diagram about the probable evolution of structural works. 


Jan 532 Jan 533 Jan 534 Jan 535 Jan 536 Jan 537 Oec 537 

site clearance 
ZZÜZD 

initial settmg-out and construction of ground level piers and colonnades 
aisle vaults and Corner ramps 

gallery level piers and colonnades to spnngs of mam arches 
gallery vaults fcorner bays probably last) 


main arches and upward contmuation of piers 

addition of projeciions to piers at ground and gallery level 
initial construction of tympana (partly rebuilt later) 


semidomes and pendentives and further contmuation of piers 


dome 

ground-level complction (mitially in the aisles oniy) 


gallery-level completion 

upper level completion and scaffold remova? 


furmshmg 


Fig. 2 - Temporal evolution of structural works according to Mainstone [ 1 ] 

The construction progress of the main structural elements of the monument at the 
early erection stages can be subdivided into three phases, which are visualized in the fol¬ 
lowing to allow a better understanding of the evolution of such a complex structural Sys¬ 
tem: 














































phase #1\ main piers are built, together with buttress piers, the linking arches and the 
piers of the secondary System; all structural elements are built from the ground level 
up to the level of the capitals of aisle’s columns. 



phase #2\ primary and secondary piers are extended up to second cornice, even if the 
transversal section is reduced with respect to the lower one. Buttress piers are 
erected too, but brick masonry is used instead of stone masonry and their building 
proceeds with increasing speed from the gallery level onward. 



phase #3\ scaffoldings are taken down from both main and secondary arches restrained 
by main piers. Horizontal thrusts Start to deform piers especially along the North- 
South direction (due to the thrust exerted by East and West arches) for mainly two 
reasons: the weakness of the link elements between main and buttress piers, and the 
differences in elastic moduli and Settlements between the two structural elements. 
Mainstone supposes that self-sustaining centerings have been used in the erection of 



minor arches, but he assigns a low probability to the fact that the same System could 
also be used for main arches; as a matter of fact, self-sustaining centerings are very 
difficult to be used when dealing with elements possessing large dimensions and, 
moreover, the presence of other structural elements at lower level could have sug- 
gested to adopt usual scaffoldings sustaining the arches during their erection. How- 
ever, this fact does not necessarily indicate that horizontal thrusts were completely 
eliminated in the construction phase (as the use of vertical scaffoldings would) as 
some reported experienced problems during the taking down of the scaffoldings 
could imply. 



Fig. 5 - Phase #3 of Hagia Sophia erection 


3.1 Deformation at the level of dome’s impost 

Some of major problems which are clearly observable at present time are pointed out 
by the large difference between diameters at the dome’s base level. With respect to the 
original circular geometry, the theoretical circumference is today deformed and diameters 
along major monument’s axes (North-South and East-West) show very different values. 

On the cornice of the present dorne two different symmetrical portions can be located 
(respectively corresponding to South and North side of the dorne) which have been built 
during the first reconstruction in 558; on these sides some reference lines carved on the 
marble cornice are present and, if one considers these curves as lying on a circular por- 
tion, the original radius of curvature can be estimated. 

Moreover, the distance between the center of the two circular sectors allows the es- 
timation of the total displacement between the top of the main piers along the North-South 
direction, due to both the flexural strain of the piers and the movements of the keystone of 
the main arches under the thrusts of the dorne itself. 

Mainstone [1] reported the distance between the two centers to a value of 2.83 m; 
some more recent analyses, after a photogrammetric survey on the monument [H], indi- 
cated in 3.338 m the present value of this distance. 

Last data Processing made by Prof. Kenichiro Hidaka and his co-workers allowed to 
demonstrate that the curves on the marble exhibit appreciable differences with respect to 
theoretical circular sectors; this fact, due to evident flexural effects on the arches, leads to 





the needing of a new estimation of the distance between centers of the two dome’s por- 
tions, once flexural effects have been subtracted from. After the new estimation procedure, 
the distance between the two semicircies can be evaluated in 2.522 -r 2.562 m [5], and this 
value can then be used as a reference vaiue in the following when deaiing with the defor- 
mation analysis of the main piers. 

Half of this value (1.271 m on average) represents the horizontal displacement of the 
keystone of each of the main arches (North and South sides) with respect to the vertical 
plane containing them, and without any contribution due to flexural effects induced by the 
final dorne in 562. 

Moreover, as it will be shown in following paragraphs when deaiing with elastic-linear 
modeling of the piers, the out-of-plane movement of the arches is roughly the double of the 
displacement of the top of the main piers. In more details, the keystone of the arch lying in 
the East-West plane (YZ plane in the following) moves in the X direction (i.e. in the vertical 
plane containing the North-South direction) by a quantity which is almost twice the move¬ 
ment of the top of the pier. 

It is then clear that the displacement between the centers of the two semicircies 
could be due to a movement of the end-section of the piers only if this movement is equal 
to half of the value of the global horizontal displacement of the keystone of the arches. As- 
suming the difference between the two centers equal to 2.542 m (average of the quantities 
suggested by Prof. Hidaka), the end section of each main pier should exhibit a movement 
equal to a quarter of this value (about 0.635 m). 

According to Mainstone’s analysis about the temporal evolution of the deformation 
(see Fig. 6), about 60% of the total deformation of the piers in the North-South direction 
should have occurred during the first collapse of the east portion of superior structures in 
558. Moreover, Mainstone observed that columns and arches placed as a strengthening of 
the main piers just after the construction of the two semidomes show an inclination that is 
roughly equal to 80 -f 85 % of the corresponding inclination of piers; is then possible to 
form the hypothesis that, in the original Situation, the inclination of piers was about 15 -f 
20% of the present one. 


Percentage of the present inclination 



Fig. 6 - Time-evolution of deformations on main piers (after Mainstone, [1]) 












By the light of previous consideration, displacements of the impost level of the arch in 
North-South direction (after the scaffoldings were taken down and betöre the construction 
of the two semidomes) were in 535 equal to about 0.11 m (15 -r 20 % of 0.635 m). 

The aim of foliowing sections is then to quantify, by means of as simple as possible 
analyses, the displacement of the end-section of piers, in order to assess whether or not 
displacements obtained in an analytical way could justify large values such as those which 
probably came through at the early phase of the monument’s erection. 

4. The “pier-system” 

According to Mainstone’s hypothesis, more representative elements of primary struc- 
tural System (i.e. piers and their buttress) exhibited their first displacements even during 
their construction. It is then necessary to focus the attention on these elements in order to 
quantify movements which could have been occurred at a very early stage of the construc¬ 
tion. In the foliowing analyses the attention will then be focused on the “pier-system”, that 
is the System composed by the main piers, their buttress and the linking elements, looking 
for the evaluation of the tendency of the System to move out of the plane of the main 
arches after having taken down the scaffoldings. 

In the foliowing models only a quarter of the whole structure has been modeled be- 
cause of the whole symmetry; erection phases are then thought to have progressed sym- 
metrically, as it has been usually assumed. A simplified scheme can be constituted by the 
“pier-system” as it has been sketched in Fig. 7. The “pier-system” is actually composed by 
two different Systems: the pier and the buttress. It is then necessary to estimate the effec- 
tiveness of the linking elements, that is the possibility that forces can be easily transferred 
between the two sub-systems. 




Fig. 7 - The “pier-system” (measures in m) 



A linear elastic analysis has been performed by using a Finite Elements Code and by 
using simple one-dimensional elements (frame elements). The structural components re- 
ported in Fig. 7 have been modeled and the analysis pointed out that the whole System 
behaves essentially as if it were composed by the internal element (pier) only; as reported 
in Table 1 (where a comparison has been made between the isolated pier and the whole 
“pier-system”), the effectiveness of the link offered by the arches is very poor, and the only 
substantial differences arise when considering the behavior along the “strong” axis (i.e. 
along direction X). In the following, the whole System can then be well represented by the 
behavior of the internal pier only. 


BASE REACTIONS 


F x 

[MN] 

F y 

[MN] 

F z 

[MN] 

M x 

[MNm] 

M y 

[MNm] 

M z 

[MNm] 

ISOLATED CANTILEVER 

1) Load -Z + s.w. 

0.00 

0.00 

13.635 

0.00 

0.00 

0.00 

2) Load -X 

1.00 

0.00 

0.00 

0.00 

24.098 

0.00 

3) Load -Y 

0.00 

1.00 

0.00 

-24.097 

0.00 

0.00 

Load combination 
D+2)+3) 

1.00 

1.00 

13.635 

-24.097 

24.098 

0.00 

“PIER-SYSTEM” 

1) Load -Z + s.w. 

-0.32 

0.00 

14.547 

0.00 

-4.296 

0.00 

2) Load -X 

0.94 

0.00 

0.00 

0.00 

22.900 

0.00 

3) Load -Y 

0.00 

0.97 

0.00 

-23.309 

0.00 

0.17 

Load combination 
1)+2)+3) 

0.63 

0.97 

14.542 

-23.309 

18.604 


% difference 

(“pier-system” vs. iso¬ 
lated cantilever) 

63% 

97% 

107% 

97% 

77% 

- 


TOP DISPLACEMENTS 


U x 

[cm] 

U Y 

[cm] 

Uz 

[cm] 

(px 

[mrad] 

Cp Y 

[mrad] 

9z 

[mrad] 

ISOLATED CANTILEVER 

Load combination 

1 )+2)+3) 

-1.53 

-1.60 

-0.122 

1.273 

1.053 

0.000 

“PIER-SYSTEM” 

Load combination 

1 )+2)+3) 

-1.29 

-1.55 

-0.131 

1.228 

0.916 

-0.022 

% difference 

(“pier-system” vs. iso¬ 
lated cantilever) 

84% 

97% 

107% 

96% 

87% 

- 


Table 1 - Comparison of base reactions and top displacements between an isolated canti- 
lever (main piers only) and the “pier-system” [applied forces are seif weight and concen- 

trated load at top with intensity equal to 1 MN] 

























































































5. Numerical modeling 

5.1 Problems related with numerical modeling 

Mainstone focused one of the possibie problems related with Hagia Sophia’s con- 
struction: the small period of time used for its construction did not allowed the mortar and 
the bricks to correctly settle, then increasing the deformation of the structure. This aspect 
can be stressed also due to the fact the Byzantine masonry is composed by mortar joints 
as thick as the bricks, then leading to an inhomogeneous mechanical behavior. 

Several measures have been adopted to control the tendency toward large deforma¬ 
tion experienced during the construction of the monument, but they somehow failed, 
probably because of the time-depending behavior of the masonry (related to creep phe- 
nomena). 

The aforementioned aspect does not allow to completeiy characterize the structural 
behavior of the monument by simply using linear elastic models; on the other hand, the 
use of complex nonlinear modeling can often iead to misleading results, mainly because of 
the large amounts of uncertain parameters to be used for a correct modeling of the non¬ 
linear mechanical and time-dependent characteristics. 

In the present study, it has been chosen not to increase the complexity of the finite 
element modeling, then leaving the numerical model as simple as possibie and within the 
limits of a linear-elastic analysis. Non linear aspects have been taken into account by 
means of simplified models of some structural elements, which have been used to evalu- 
ate linear-elastic equivalent properties to be used in the FE models. In following Para¬ 
graph, some details will be given about the mechanical characterization of the materials 
adopted in the subsequent analyses. 

5.2 Materials parameters 

Piers cross section is varying along the height of the pier itself; in the lower part (up 
to a level of about 18 m) the masonry is made by stones and a filling with cohesive proper¬ 
ties, while in the upper part the masonry is simply made by mortar and bricks. 

Transversal cross section of piers has been preliminarily brought back to an equiva¬ 
lent rectangular section, as indicated in Fig. 8a, while data derived from literature have 
been used to characterize the mechanical properties of the masonry. 

• Stone masonry (lower levels) 

The thickness of the external layer made by stone masonry is not well known, as well 
as their mechanical properties and those of the filling material. In the view of looking 
for the values of the parameters describing the mechanical behavior of an “equivalent” 
homogenous material, following assumptions have been made. 

As a first approximation, once values for Young’s moduli and areas have been chosen 
for the two materials ( E-i for stone masonry whose area is equal to A 1; E 2 for the filling 
with area A 2 ) the “equivalent” Young’s modulus can be chosen by imposing vertical 
displacement equal for the two materials, once a load has been applied. Indicating 
with N the global vertical force and with Ni and N 2 the amount of the force on the stone 
masonry and the filling respectively, the equality of vertical strains implies that 

N N 1 _ N 2 

E eq ' t A 1 + A 2] E 1 A 1 E 2 A 2 
then obtaining 



NE 1 ■ A 1 


e N 


N E 2 a 2 


Ni = 


E-) ■ A-| + E 9 ■ Ac 


2 - 


E-| ■ A-| + E 9 ■ Ac 


The denominator represents the whole stiffness of the cross-section, so that the 
equivalent modulus can be estimated as 


E-| ■ A-| + E 2 ■ A 2 

eq= 



Fig. 9 - Transversal cross-section of the pier (measures in m): 
a) geometrical characteristics; b) layers constituting the piers in the lower part 


E eq value is then dependent on the ratio between the two areas, the latter depending on 
the thickness of the external stone layer. By assuming different values for this layer, 
values reported in Table 2 can be obtained, where E^IOOOO MPa and E 2 =2000 MPa 
have been adopted (according to next Table 3). 

Because of the lack of correct information about the effective thickness of the stone 
layer, a value of E eq =5000 MPa has been chosen in the following, close to the obtained 
average value. 


Stone thickness 

Eeq 

[m] 

[MPa] 

0,3 

3488 

0,6 

4832 

0,9 

6032 

average 

4784 


Table 2 - E eq values as a function of the stone thickness 













Author 

Stone 

masonry 

Mortar 

Building 

Numerical 

modeling 

Croci et al. [6] 

10000 

2000 

Flagia Sophia 
(Istanbul) 

Flidaka et al. [7] 

10000 

10000 

Flagia Sophia 
(Istanbul) 

Laboratory 

tests 

Penelis, 

Karaveziroglou [8] 

- 

2450 

Rotunda of Thessa¬ 
loniki 

(early Byzantine age) 

Theocharidou [9] 

- 

2000 

Flagia Sophia in 
Thessaloniki 
(beginning of VII 
Century A.D.) 


Table 3 - Young’s modulus values (in MPa) for different masonry 

(from cited bibliography) 


Stone masonry (upper levels) 

As a first approach, mechanical characteristics of the masonry composed by alternate 
mortar and bricks layers with equal thickness could be determined by considering a 
pure compressive stress state and neglecting all the effects related to strains orthogo¬ 
nal to the direction of the applied force. By utilizing a scheme in which the two different 
layers can be represented by a series of springs with different characteristics which, 
under a load P, exhibit a displacement proportional to their stiffness ki e k 2 (Fig. 9). 
Being the stiffness of bricks Iayer higher than the stiffness of mortar layers (ki » k 2 ), 
the resulting equivalent stiffness is then 

K ~ k 2 

As a consequence, it can be thought that the whole stiffness is governed by the 
weaker element (the mortar), especially in a first phase when the hardening of the 
mortar took place. 

According to results reported by recent researches made by the Politecnico di Milano 
[10], the Young’s modulus for the brick masonry can then be assumed equal to the 
lower value reported in Table 3, that is E m = 2000 MPa. 



Fig. 9 - Series of springs: equivalent stiffness K 




















5.3 Linear modeling 

The characteristics of the materials singled out in previous paragraph can be used for 
the determination of the stress-strain state in a liner-elastic condition; the analyzed struc- 
ture is composed by a quarter of the whole structure before the erection of the dome, that 
is by the “pier-system” and by two semi-arches in the North-South plane (XZ) and in the 
East-West plane (YZ) respectively. The restraint offered by the external (non modeled) 
zones has been taken into account by inserting appropriate symmetry restraints. 

The portion of the structure has been modeled by using “beam type” elements (see 
Fig. 10); constituting materials possess the characteristics reported in Table 4. Obtained 
displacements and internal forces under vertical loads only (seif weight) are reported (for 
the most representative points) in next Table 5. 


Parameter 

Value 

Equivalent Young’s modulus for stone masonry 

E s = 5000 MPa 

Young’s modulus for brick masonry 

E b = 2000 MPa 

Poisson coefficient (both materials) 

v=0.3 

Specific weight (both materials) 

y=18 kN/m 3 


Table 4 - Mechanical characteristics used in the linear elastic models 


pos (B): keystone 
of (YZ) arch 




Fig. 10 - Linear elastic model: deformed configuration 


5.4 Non linear modeling 

In previous paragraph it has been pointed out that the “pier-system” can be ade- 
quately schematized by a cantilever constituted by the internal pier only; non-linear calcu- 
lus can then be carried out for an isolated cantilever element, looking for an extrapolation 
of information to be used in subsequent analyses. 























In following analysis a numerical nonlinear analysis method has been used, based on 
the one which has been originaliy set up at the Civil Engineering Department of the Uni- 
versity of Florence (Angotti et al. [11]). 

By means of the numerical model (see Fig. 11), it is possible to determine the stress- 
strain state of a slender element under axial and flexural loads, by taking into account dif¬ 
ferent sources of nonlinearities, both geometrical (linked to P-S effect) and mechanical 
(nonlinear stress-strain relationships); influence of shear deformation can be taken into ac¬ 
count, too. 


Displacements 

U x 

[cm] 

U Y 

[cm] 

Uz 

[cm] 

<px 

[mrad] 

Cp Y 

[mrad] 

9z 

[mrad] 

pos. (A): impost level 
of arches 

-1.58 

-1.57 

-0.25 

1.180 

-1.111 

-0.033 

pos. (B): keystone of 
(YZ) arch 

-3.31 

0.00 

-3.69 

0.000 

-1.111 

-0.033 

pos. (C): keystone of 
(XZ) arch 

0.00 

-3.51 

-0.48 

1.180 

0.000 

-0.033 

Internal forces 

F z 

[MN] 

Fz 

[MN] 

F z 

[MN] 

M x 

[MNm] 

M y 

[MNm] 

M z 

[MNm] 

pos. (A): top level of 
internal pier 

-5.060 

-0.999 

-0.793 

0.737 

-1.700 

0.000 


Table 5 - Results from linear-elastic model (remark: the out-of-plane displacements of 
keystone points of arches is about the double of the displacement at impost level) 



MAIN HYPOTHESES 

• Variable section cantilever, with 
non uniform mechanical proper- 
ties 

• Loads applied at the top level 
section (in every possible direc- 
tion) and loading due to seif 
weight 

• Finite element discretization of 
the whole structure 

• Cross-section discretization by 
means of small rectangular ele- 
ments 

• Geometrical nonlinear behavior 
(P-S effect) 

• Mechanical nonlinear behavior 
(material with no tensile stresses) 

• Shear deformation accounted 


Fig. 11 - Nonlinear analysis: reference model and main hypotheses 
































The adopted makes use of the equations reported in Tabie 6, where the cantiiever 
with height L is subdivided into n finite elements of equal height Az=L/n. The solution pro- 
cedure is an iterative one: the analysis is stopped when convergence is reached between 
the stress and the strain state in all the elements. 


1) RELATIONSHIP BETWEEN DISPLACEMENTS AND CURVATURES 


v i 


Az 2 1 


2 r 0 
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v,: horizontal displacement 
of i-th finite element 
(numbering Starts from 
the clamped end) 

n: radius of curvature of i- 
th finite element 

r 0 : radius of curvature of 0- 
th (base) finite element 

Av iit : horizontal displace¬ 
ment of i-th finite ele¬ 
ment due to shear de- 
formability 

z\\ level of i-th finite element 


2) RELATIONSHIPS BETWEEN BENDING MOMENTS AND DEFORMED 
CONFIGURATION 


M n = M 

M n -i =M + H-(L-z n _ 1 )+P-(v n -Vn-i) 

APi: axial force on i-th finite 


element 

n-1 

Mi: bending moment at 

M i =M + H(L-z i )+P(v n -v i )+ £AP k .(v k - Vi ) 

lower level of i-th finite 

k=i+l 

element 

n 

M 0 : base bending moment 

Mq = M + H-L + P-v n + AP^ v k 

(clamped end reaction) 

k=l 



3) GEOMETRICAL CHARACTERISTICS OF RESISTANT SECTION 


Ac= . 

j' dA 

Sx=J 

y dA 

For each finite element the 
Integration domain extends 

s y =j 

i 

xdA 

9 „ 

^xy — 

f 

x ■ y dA 

to rectangular elements with 
negative (compressive) 

J x=j 

y~ dA 

J y =j 

x 2 dA 

stress 


4) RELATIONSHIP BETWEEN INTERNAL FORCES AI 

ND CURVATURES 
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L J J 
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K x : curvature in (YZ) plane 
K y : curvature in (XZ) plane 
£z: axial deformation 


Tabie 6 - Main relationships used in the adopted nonlinear calculation 






Deformed configuration is defined by means of the centroid displacements in the two 
planes (XZ) and (YZ), being X and Y the directions corresponding to main inertia axes of 
the cross-section. Starting from undeformed geometry, the procedure Starts by evaluating 
the values of the internal forces (ensuring global equilibrium through eqn. #2) and then of 
internal stresses (eqn. #3). New values for internal forces are then evaluated by using the 
stress-strain relationships (eqn. #4) and hence the displacements (eqn. #1), in both main 
structural planes. The procedure stops when the obtained stress state and the values of 
internal forces are in accordance with the strain state and displacement of the centroid of 
the cross-section. 

In the algorithm it is possible to decide whether or not the different sources of nonlin- 
earities have to be taken into account, so estimating the mutual influence among all the 
different factors involved in the evaluation of the global displacement. 

Elastic moduli are the same adopted in the linear-elastic case while for transversal 
modulus of elasticity G a different value has been used. In order to take into account the 
high shear deformability due to the thickness of mortar layers, a value equal to G=E/5 has 
been assumed, then abandoning the linear-elastic relationship. 

Loads applied at the end section of the cantilever beam are those derived from lin¬ 
ear-elastic calculations performed over a quarter of the whole structure (reported in the 
last row of previous Table 5). 

In the following, main obtained results are reported; Table 7 reports the displacement 
values at the top level obtained by three different analyses: by considering the no tensile 
stress condition only, by considering both no tensile stress condition and P-S effect, and by 
considering the previous two nonlinear effects plus the shear deformability. 


Performed calculation 

U x 

[cm] 

U Y 

[cm] 

V(Ux 2 +U y 2 ) 

[cm] 

L) Linear elastic 

-1.71 

-1.60 

2.34 

1) No tensile stress condition 

-2.85 

-2.41 

3.73 

2 ) No tensile stress condition + 
P-S 

-3.47 

-2.27 

4.15 

3) No tensile stress condition + 
P-S + shear deformability 

-3.71 

-2.41 

4.42 


Table 7 - Displacement at cantilever top level 
(E s =5000 MPa, E b =2000 MPa) 


From reported results it can be noticed that the nonlinear effect giving the highest 
contribution is the one relate to nonlinear mechanical behavior (absence of tensile resis- 
tance); by inserting this effect, displacements increase because of the increased curvature 
of each single finite element. 

The ratio between the top level displacement obtained by the performed nonlinear 
analysis and the same value when dealing with a linear-elastic calculation is about 2.0. 
This fact implies that, if one would have only linear-elastic analyses to be performed, dif¬ 
ferent elastic properties of the materials should be adopted with respect to the initial one in 
order to take (indirectly) into account nonlinear effects. “Correct” moduli values should be 
as low as 50% of the initial ones, that is equal to 




• stone external masonry + filling: E s = 2500 MPa 

• brick masonry: E b = 1000 MPa 

The obtained values for elastic moduli can be used for modeling the elastic behavior 
of the “pier-system” inserted in the linear elastic model of one quarter of the whole struc- 
ture, leaving values of elastic moduli of the other structural elements (arches) at the same 
values as before. linder the same vertical actions, internal forces at the top level of the 
pier-system are different with respect to those reported in Table 5; by introducing the new 
obtained top-level forces into the nonlinear procedure sketched before, the new obtained 
displacement values for the top-level section would be as reported in Table 8. 


Performed calculation 

U x 

[cm] 

U Y 

[cm] 

V(u x 2 + U Y 2 ) 

[cm] 

1) No tensile stress condition 

-1.50 

-1.33 

2.01 

2 ) No tensile stress condition + 

P-S 

-2.18 

-2.04 

2.98 

3) No tensile stress condition + 

P-S + shear deformability 

-2.61 

-2.02 

3.30 


Table 8 - Displacement at cantilever top level 
(E s =2500 MPa, E b =1000 MPa) 


The ratio between the displacement of the top-level of cantilever obtained by a non¬ 
linear analysis and the same quantity obtained by a liner-elastic calculation is now equal to 
1.7, so that “correct” values for elastic moduli should be 

• stone external masonry + filling: E s = 2941 MPa 

• brick masonry: E b = 1177 MPa 

The procedure should now be repeated until convergence; however, stating the fact 
the obtained values are intermediate with respect to those initially assumed and those ob¬ 
tained after the first Step of this procedure, one would expect that new analyses would lead 
to top-level displacement values included in the ränge delimited by previous values, that is 
to an amplification factor (with respect to linear-elastic modeling) within the ränge 1.7^2.0. 

As it has been noticed in previous paragraphs, displacement along North-South di- 
rection at the impost level of main arches placed in East and West sides of the building 
should be, after the erection of the two semidomes, as high as about 11 cm in order to jus- 
tify the measured difference in the main diameters of the dorne. 

In the reported analyses, performed by using an iterative linear-nonlinear simplified 
approach, the maximum obtained displacement for the same section (in an erection stage 
corresponding to the Situation immediately before the erection of the two semidomes) can 
be assumed as within the ränge 3^4.5 cm. 

6. Concluding remarks 

In the present work some initial evaluation of deformation happened in some struc¬ 
tural elements during early erection stages of Hagia Sophia have been reported. In more 
details, the displacements at impost level of the dorne have been analyzed by investigating 
the behavior of main piers, in order to explain the initial “ovalling” tendency in the circular 
cross-section of the dorne. 




Calculations have been performed by means of very simple linear models and by it- 
erating a nonlinear procedure based on the behavior of main elements only; the nonlinear 
analysis shown that simple linear-elastic models could lead to evident underestimation of 
displacements which have probably occurred in main piers. 

The obtained results, however, did not seem to completely justify the actual values 
which should have occurred in a very early stage of the monument erection. Further 
analyses will then be performed in order to take into account also viscous phenomena, so 
to include an effect which probably could lead to more realistic displacements’ values, 
closer to the ones which took place in Hagia Sophia before the erection of the main dorne. 
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